Coulomb drag by small momentum transfer between quantum wires 
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We demonstrate that in a wide range of temperatures Coulomb drag between two weakly coupled 
quantum wires is dominated by processes with a small interwire momentum transfer. Such processes, 
not accounted for in the conventional Luttinger liquid theory, cause drag only because the electron 
dispersion relation is not linear. The corresponding contribution to the drag resistance scales with 
temperature as T 2 if the wires are identical, and as T 5 if the wires are different. 
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Electrons moving in a conductor generate a fluctuat- 
ing electric field around it. This field gives rise to an un- 
usual transport phenomenon. Coulomb drag between two 
closely situated conductors [l| . The structure of the fluc- 
tuating field is determined by electron correlations within 
the conductors. Correlations are stronger in conductors 
of lower dimensionality. Tomonaga-Luttinger model cap- 
tures some aspects of the correlations in the case of one- 
dimensional conductors (quantum wires). Within this 
model, the Coulomb drag was studied in 2, 3]. In a typi- 
cal setup jj, |a] a dc current flows through the active wire 
1, while the bias applied to the passive wire sets 1% = 0, 
see Fig. ^ The drag resistivity (drag resistance per unit 
length of the interacting region) is then denned as 

r = — lim — — 77 — — r— — . (1) 
Zi— o 2whL dh 

The only source of drag in the Luttinger liquid is in- 
terwire backscattering, associated with a large momen- 
tum transfer between the wires. The model predicts a 
distinctive temperature dependence of the corresponding 
contribution r2k F to the drag resistivity Q). In the case 
of identical wires r 2 k F oc l2k F eA ^ T at the lowest temper- 
atures 0, . Here l2k F is the scattering length charac- 
terizing the interwire backscattering. At temperatures 
T above the gap A, this exponential dependence is re- 
placed by a power-law, r2k F ~ l^k (T/tF) 1 ~ 7 , where cf 
is the Fermi energy. The exponential temperature depen- 
dence of r2k F indicates a formation of a zig-zag charge 
order due to the 2/cF-component of the interwire interac- 
tion |2|, yj ■ To the contrary, the exponent 7 > in the 
power-law portion of the function r 2 k F (T) is determined 
by the interactions within the wires; 7 = in the absence 
of interactions |(| ■ This renormalization of T2k F is similar 
in origin to the suppression of the conductance of a Lut- 
tinger liquid with an impurity : in both cases repulsive 
interactions enhance the backscattering probability when 
temperature is lowered. 

However, forward scattering between the wires also 
induces drag. To see this, one has to go beyond the 
Tomonaga-Luttinger model and account for the nonlin- 
earity of the electronic dispersion relation. If the electron 



velocity depends on momentum, then even small (com- 
pared to 2kp) momentum transfer results in drag. 




FIG. 1: Coulomb drag between quantum wires. A dc current 
1 1 flows through the active wire (1). A voltage bias V2 is 
applied to the passive wire (2) in such a way that I2 = 0. 

The small (ro) and large (r2k F ) momentum transfer 
contributions to the drag are inversely proportional to 
the corresponding scattering lengths Iq and l2k F , respec- 
tively. Their ratio lo/hk F cx e~ ikpd depends strongly on 
the distance d between the wires. If the drag effect is 
employed to study the correlations within a wire, rather 
than the zig-zag order induced by interwire interaction, 
then d must be large: kpd > 1. In this case the gap 
A ~ £F(fc_F^2fc F ) -1 ^ 7 becomes narrow, and the role of 
forward scattering increases. 

In this paper we demonstrate that the drag resistivity 
between weakly coupled wires is dominated by the for- 
ward scattering in a wide temperature range. Even for 
identical wires, which is the most favorable for backscat- 
tering case, ro oc T 2 wins over r2k F at all T above 
T* - € F {lo/l 2 k F ) 1/{1+l) - For different wires, T2k F is ex- 
ponentially small at T < u5n, whereas ro has a power-law 
low-temperature asymptotics; here Sn is the mismatch of 
the electron densities between the wires and u is the char- 
acteristic plasma velocity; hereafter we set K = 1. 

The Hamiltonian of the system, H = H\ + H2+ H12, is 
given by the sum of the Hamiltonians of the two isolated 
wires i = 1, 2, 



II, = I dxi/;j(x)-^ipi(x) + Hint, p=-id/dx, 



H int = / dxdypi(x)Uu(x-y)pi(y), 



(2) 



2 



and of the Hamiltonian of the interwire interaction, 



H\2 



dx X dx 2 p\(xi)Ui 2 (xi -X 2 )p 2 (x 2 ); (3) 



here pi(x) = ipl(x)i/ii(x). 

We will assume that the interactions are screened by 
the nearby metallic gates. The screening length oL is 
set by the distance to the gates and is typically 0, Q of 
the order of the separation d between the wires, d s ~ d. 
The short-distance cutoff dy of the Coulomb potential 
is provided by the wire width do for i — j or by the 
interwire distance d for i 5^ j. The Fourier transforms 
Uij(k) = J dxe tkx Uij(x) are rapidly decreasing func- 
tions of k with U io {k) oc e~^ d " for \k\ » l/d s , and 
Uij(k) rs const for \k\ <gC Note that dy ~ d s ~ d 

for interwire interaction Ui 2 (k). Thus, its fc-dependence 
is characterized by a single scale fco ~ l j d. 

Because of the interaction H X2 , electrons in the wire 2 
experience a force |2j whose density is given by 



T 2 = \ dx{dU 12 {x)/dx){ Pl {x)p 2 {Q)). 



(4) 



Since there is no current in the wire 2, this force must 
be balanced by an external electric field, en 2 £ 2 +^2 = 0, 
where £ 2 — V 2 /L and = (pi) is the concentration of 
electrons in the wire i. At T > A (see the discussion 
above) the correlation function in the r.h.s. of Eq. (0J 
can be evaluated in the first order in Ui 2 , 



T = —\w% *E&(*)Si(*,«)&(-*,-«), (5) 
L en 2 J (2uy 



where Si(k 1 u>) are the dynamic structure factors, 
Si(fc,w) = fdxdt e iut ~ ikx ( Pi (x,t) Pi (0,0)) , 



calculated in the presence of a finite current I\ in the 
wire 1. The structure factor S2(k,u>) in the wire 2 coin- 
cides with its equilibrium value, S 2 (k,u>). The electronic 
subsystem in the wire 1 is in equilibrium in the refer- 
ence frame moving with the drift velocity Vd = I\/en\ 
in the direction of the current. Therefore the structure 
factor S\ is obtained from the equilibrium value Si using 
the Galilean transformation: Si(k, ui) — Si(k,ui — qvd)- 
Equations and © then yield 



= / dkduj 



k 2 U 2 2 (k) dSi{k,Lu) 



S 2 (-k,-cj). (6) 



87r 3 niri2 duj 
Now we use the fluctuation-dissipation theorem, 

2Ai(A,w) 



Si(k,uj) = 
to further simplify Eq. (|^|). 



1 - e-"/ T 



rdk Td, y U (*M-(*;^M, (7) 



'0 jo 



4ir 3 nm 2 T sinh 2 (w/2T) 



Here Ai is the imaginary part of the retarded density- 
density correlation function; Ai(k,uj) = Ai(—k,uj) = 
— Ai(k, — cj). Equation JJJ) was derived by different 
means in |£| ; similar expressions have been also obtained 
for nonintcracting systems with disorder |IJ. Here we 
demonstrated the validity of Eq. Q for clean interacting 
systems. 

We start with the evaluation of the drag resistivity for 
nonintcracting electrons (Uu = 0). Concentrating on the 
small momentum transfer contribution to r, we consider 
the limit lo/l 2 k F — > 0, thus setting T2k F = 0. In this case 
the main contribution to the integral over k in Eq. Q 
comes from small momenta k kp and small energies 
uj <C £f ■ At these values of k and ui, functions Ai(k,uj) 
are sharply peaked at u>i = Vik, where Vi = irrii/m are 
the Fermi velocities in the two wires. For a given k < 2fcp 
the widths of the peaks can be estimated as 



Su>(k, T) — max |fc 2 /m, kT/kp} 
Equation (jHJ and the exact f-sum rule, 



f , . . j . irn, k 2 

/ duj ujAi(k,uj) = — , 

Jo m 2 



(8) 



(9) 



allow us to estimate the peak heights: Ai ~ k/2Su>. If 
the difference between the Fermi velocities is small, 



Sv = \vi — v 2 \ 
then Eq. reduces to 
1 



Ttn/m, n = (ni + n 2 )/2, 



°° dk —^M^— a (k,T), (10) 



sinli (v F /c/2T) 



a(k,T)= diuAi(k,uj)A2(k,Lo). 
Jo 



(11) 



The function a(k,T) depends on 5v. If the wires are 
identical (Sv = 0) , then Eq. pi(l and the above estimates 
for Ai yield 



a(k, T) 



4Su(k,T)' 



(12) 



There are two competing scales in the integrand of 
Eq. I|10[) . The first scale, fco ~ 1/d <C kp, character- 
izes the fc-dependence of the interwire interaction Ui 2 (k). 
The typical wave vector of thermally excited electron- 
hole pairs, T/vf, defines the second scale. The two scales 
coincide at T — Tq — vpk^. At T <C To one can replace 
U\ 2 (k) by t/i2(0) in Eq. (fTTj|) . Furthemore, we use a in 
the form of Eq. 1)12(1 at T — 0, which results in 



£1 
h 
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1 


~u 12 (oy 




To ~ 


2ttvf 



(13) 



withci = 7r 4 /12. Use of exact form of Ai(k, u>) inEq. I|ll|) 
changes only the numerical coefficient, ci = 7r 2 /4. 



The increase of temperature T above To results in a 
saturation of the drag resistivity. Indeed, at T C T <C 
e F one can expand sinh(upfc/2T) in Eq. (|1U|I and use 
Slu = kT/k F for the peak width in lj*H21) . This yields 



kdk U 2 2 (k) 



1 

T 



(14) 



Further increase of T leads to the decay of the drag, 

rcxl^iTo/epfiT/ep)- 3 / 2 , T » e F , (15) 

similar to the two-dimensional case [To| . 

We now consider wires with slightly different Fermi ve- 
locities 5v > 0. In this case the peaks of Ai(k, uj) are sep- 
arated in u> by k8v. We define a new temperature scale 
Ti = kp5v by equating the separation to the peak width 
©■ We assume this scale is small, Ti <C T . The differ- 
ence between velocities does not affect the drag at T ^ 
T\. However, at T <C Xi the drag resistivity is suppressed 
exponentially. To obtain the leading asvmptotics of r(T) 
it is sufficient to use the T = limit 11] of Ai(k,ui) in 
Eq. (HTJ, a(k, 0) = (m/Ak) (k - m5v) 9(k - mSv). Equa- 
tion i|l(J|) then results in 



(p 



(16) 



The activational temperature dependence Eq. l|ltj|) 
holds for all T <C Ti, because for noninteracting elec- 
trons at T = the product A\{k, uj)A 2 (k, w) is exactly 
zero ^lj at k < mSv. If electrons interact, some overlap 
of A\ and A 2 exists even at small k <C mSv. This yields 
a further contribution to r, that has a power-law temper- 
ature dependence. We will evaluate this contribution for 
weak intrawire interaction. 

It is convenient to write A(k,u>) (we suppress the 
index i in the following) in the form A{k,u>) — 
[S(k,ui) — S(—k, — cj)] /2 and use the Lehmann represen- 
tation for the dynamic structure factor: 

S(k,u)= ^-Y.^ n \P^)\ 2 ^-E n + E gs ). (17) 

n 

Here L is the system size, \gs) is the ground state, and 
Pk = SpV'p+feV'p- We evaluate the matrix element in 
Eq. ((1711 in the first order in the intrawire interaction. 
The nonvanishing at to — Vpk 3> 8u> contribution re- 
sults from the processes in which the unperturbed fi- 
nal state \n) in Eq. i|17|) has two electron- hole pairs: 
In}' -* = i^p+qi^pt/jpi _g<Vy |0). This contribution is 



SS(k,oj) — — [ dpdp'dqdq' S(q — q' — k) 
^ J 



(18) 



X 6(W - £, p+q +£,p~ £,p'-q> + £p') 

x /p(1 - fp+q)fp>(l ~ f P '-q')K 2 (p,p',q,q\uj), 



where f p are the Fermi functions, £ p = p 2 /2m, and 

U(q')-U(p-p' + q') U{q')-U{p-p' + q) 



K 



u — £.p+q + £p+<?' 

+ [p^p',q^ -q']. 



U — £,p+q-q' + £p 



Note that Eq. I|18|). unlike Eq. (|17|l . accounts for a finite 
temperature. At u> <C ep and fc <§; /cp, Eq. I|18|) yields 
the interaction-induced correction to A(k, u), 



8A{k,w) 



U 2 fc 4 9(lu - v F k) 
vp m 2 ui 2 — vpk 2 ' 



(19) 



where U = [U(0) - U{2k F )}/2irv F < 1. This re- 
sult is valid for u> <C ep, k <C hp, \u> — vpk\ 3> 
max{[/t> F fc, (5w(fc, T)}, and describes w) outside the 
interval Eq. (JSJ. The limit of linear electron dispersion 
relation (m — > oo) corresponds 01 to <L4(fc, = 0. 

We use Eq. (|19fl to evaluate the interaction-induced 
correction Sr to the drag resisitivity between non- 
identical wires with T\ — k F 5v e F U . At the lowest 
temperatures, Eqs. 1)10(1 and (| 1 1 1> yield 



In 



e F 



(20) 



With increasing temperature, the r(T) dependence 
changes from Eq. I|20|) to the activation law 1(16(1 . At 
T T\ the difference between wires does not affect r(T). 

We will argue now that intrawire interactions do not 
change the quadratic temperature dependence of r(T) 
at Ti « T < T , see Eq. (JTSJ. At these temper- 
atures, an estimate equivalent to Eq. 1)13(1 reads r ~ 
\U 2 2 (0)\vp 3 5uj(k T ,T) and yields r cx T 2 ; here k T ~ 
T/vp <C fci? is the wave vector of a typical electron- 
hole excitation. Interaction apparently does not affect 
the functional dependence of 8u> on k and T; the esti- 
mate |[5J| still can be used, although the coefficients 1 /m 
and 1//cf in it are affected by the interaction. 

The Tomonaga-Luttinger model is insufficient for the 
evaluation of Slu in the presence of interaction: it implies 
linear electron spectrum, which yields \l2\ 8ui = 0. Ac- 
counting for the curvature of the electron spectrum com- 
plicates the treatment of the interaction greatly. The 
width Slu can be explicitly evaluated in the Calogero- 
Suthcrland model which is characterized by a very spe- 
cific interaction potential, 



Uu(x) 



2tt 2 A(A - 1) 



mL 2 



sin nx 



(21) 



The parameter A here is related to the conventional in- 
teraction parameter g of the Luttinger liquid: g = 1/A. 
This relation follows from the definition g = v F /u in 
terms of the velocity of the collective mode (plasmon) u, 
and its value u = (jm/m)\ in the Calogero-Sutherland 
model [T3L ITil ]. For the rational values of A and at 



4 



T = the density-density correlation function is known 
exactly [l3l Hif . Due to the integrability of the model, 
Ajf fc, u>) ^ only in a finite interval of lo around lo = 
Uik |l5j. We found this interval for k < 2imi\ 



k 2 k 2 
(Vff)^— < u-u, t k < — , 



which yields for the width 



5w(k,0) 



1+gk 2 
2g m 



(22) 



(23) 



In order to estimate r we note that Eq. (J7J and the 
sum rule © remain valid in the presence of interactions 
within the wires. This allows us to follow the steps that 
led to Eq. (|13|) . Replacing vp by the plasma velocity u 
in Eq. (fT7Tf> and using Eq. f2l5|) , we find 



9 



1 + 3' 



(24) 



which agrees with our expectation for the r(T) depen- 
dence. We are not aware of a reliabale theory of Ai(k, to) 
beyond the exactly solvable case. However, the self- 
consistent Born approximation results |16| allow us to 
corroborate the estimate Soj oc k 2 jra for the peak width, 
so, apparently, the r oc T 2 dependence is universal. 



rp-3/2 




FIG. 2: Sketch of the temperature dependence of the drag 
resistivity between identical wires. The small momentum 
transfer contribution considered in this paper dominates at 
T > T*; the ratio T* /ep is exponentially small for kpd > 1 . 

First observations of drag between ballistic quantum 
wires appeared recently 0, Q • In a limited temperature 
interval, 0.2 K < T < 0.9 K, a three-fold drop in the drag 
resistance was observed Q with the increase of temper- 
ature. This drop was fit to a power-law r oc T~ 77 and 
interpeted as evidence of the Luttinger liquid behavior. 
However, the Fermi wave vector in the wires of Ref. |4| 
was estimated to be kp — 6 x 10 4 cm -1 , which yields 
€ F = h 2 k 2 F /2m* « 0.2 K (we used here m* = 0.068 m 
known for GaAs). It thus appears that the measurements 
of Ref. correspond to a non-degenerate or weakly de- 
generate regime incompatible with the Luttinger liquid 



description. An alternative explanation of the observa- 
tions 0, 0] is provided by our theory. Indeed, using the 
values of kp and d — 200 nm of we find kpd = 1.2. 
Under this condition the small momentum transfer con- 
tribution dominates at T > T* , see Fig. |2 The ob- 
served H behavior of r(T) may correspond to the 
crossover regime between the limits r(T) = const and 
r(T) oc T~ 3 / 2 presented by Eqs. JUj) and fT5|) . 

To conclude, the small momentum transfer contribu- 
tion dominates Coulomb drag at almost all temperatures 
if the distance between the wires exceeds the Fermi wave- 
length, see Fig. [21 Drag by small momentum transfer is 
possible because electron dispersion relation is not linear, 
and therefore can not be accounted for in the conven- 
tional Tomonaga-Luttinger model. 
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